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Abstract. We investigate the quantum phase transition of itinerant ferromagnets. It is 
shown that correlation effects in the underlying itinerant electron system lead to singularities 
in the order parameter field theory that result in an effective long-range interaction between 
the spin fluctuations. This interaction turns out to be generically antiferromagnetic for 
clean systems. In disordered systems analogous correlation effects lead to even stronger 
singularities. The resulting long-range interaction is, however, generically ferromagnetic. 

We discuss two possibilities for the ferromagnetic quantum phase transition. In clean 
systems, the transition is generically of first order, as is experimentally observed in MnSi. 
However, under certain conditions the transition may be continuous with non-mean field 
critical behavior. In disordered systems, one finds a very rich phase diagram showing first 
order and continuous phase transitions and several multicritical points. 
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1 Introduction 

Quantum phase transitions are phase transitions that occur at zero temperature 
as a function of some non-thermal control parameter. The fluctuations that drive 
these transition are of quantum nature rather than thermal in origin. Among the 
transitions that have been investigated arc various metal-insulator transitions, the 
superconductor-insulator transition in thin metal films, and a variety of magnetic 
phase transitions. Quantum phase transitions have attracted considerable attention 
in recent years, in particular since they are believed to be at the heart of some of 
the most exciting discoveries in modern condensed matter physics, such as the local- 
ization problem, various magnetic phenomena, the quantum Hall effects, and high- 
temperature superconductivity 

One of the most obvious examples of a quantum phase transition is the transition 
from a paramagnetic to a ferromagnetic metal that occurs as a function of the exchange 
coupling between the electron spins. The experimentally best studied example of such 
a transition is probably provided by the pressure-tuned transition in MnSi MnSi 
belongs to the class of so-called nearly or weakly ferromagnetic materials. This group 
of metals, consisting of transition metals and their compounds such as ZrZn2, TiBe2, 
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Fig. 1 Phase diagram of MnSi. The insets show the behavior of the susceptibility close to 
the transition. (After j^). 



NiaAl, and YC02 in addition to MnSi are characterized by strongly enhanced spin 
fluctuations. Thus, their ground state is close to a ferromagnetic instability which 
makes them good candidates for actually reaching the ferromagnetic quantum phase 
transition experimentally by changing the chemical composition or applying pressure. 

At ambient pressure MnSi is paramagnetic for temperatures larger than Tc — 30 K. 
Below Tc it orders magnetically. The order is, however, not exactly ferromagnetic but 
a long-wavelength (190 A) helical spin spiral along the (111) direction of the crystal. 
The ordering wavelength depends only weakly on the temperature, but a homogeneous 
magnetic field of about 0.6 T suppresses the spiral and leads to ferromagnetic order. 
One of the most remarkable findings about the magnetic phase transition in MnSi is 
that it changes from continuous to first order with decreasing temperature as is shown 
in Fig. ^ Specifically, in an experiment carried out at low pressure (corresponding to 
a comparatively high transition temperature) the susceptibility shows a pronounced 
maximum at the transition, reminiscent of the singularity expected from a continuous 
phase transition. In contrast, in an experiment at a pressure very close to (but still 
smaller than) the critical pressure the susceptibility does not show any sign of a di- 
vergence at the phase transition. Instead, it displays a finite discontinuity suggestive 
of a first-order phase transition. 

A related set of experiments is devoted to a phenomenon called itinerant elec- 
tron metamagnetism. Here a high magnetic field is applied to a nearly ferromagnetic 
material such as Co(Sei_2;Sa;)2 or Y(Coi_2;Al2:)2 At a certain field strength the 
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magnetization of the sample shows a pronounced jump. This can easily be explained 
if we assume that the free energy as a function of the magnetization has the triple- 
well structure characteristic of the vicinity of a first-order phase transition. In zero 
field the side minima must have a larger free energy than the center minimum (since 
the material is paramagnetic in zero field). The magnetic field essentially just "tilts" 
the free energy function. If one of the side minima becomes lower than the center 
(paramagnetic) one, the magnetization jumps. 

In the literature the first-order transition in MnSi at low temperatures as well 
as the itinerant electron metamagnetism have been attributed to sharp structures in 
the electronic density of states close to the Fermi energy which stem from the band 
structure of the particular material. These structures in the density of states can lead 
to a negative quartic coefficient in a magnetic Landau theory and thus to the above 
mentioned triple-well structure. 

In this paper it will be shown, however, that the two phenomena are generic since 
they are rooted in the universal many-body physics underlying the transition. There- 
fore, they are predicted to occur for all nearly or weakly ferromagnetic materials irre- 
spective of special structures in the density of states. In this paper we will emphasize 
the basic physics behind these findings, a technical derivation can be found elsewhere 

The paper is organized as follows: In Sec. |2|we sketch the derivation of an order 
parameter field theory for the ferromagnetic quantum phase transition of itinerant 
electrons, starting from a microscopic description of an interacting electron system and 
analyze its properties. In Sec. ^ we compare the possible scenarios for the quantum 
phase transition in detail, while Sec. ||is devoted to the influence of quenched disorder. 
We conclude in Sec. ^. 

2 Order parameter field theory 

In a pioneering paper that was the first application of the modern theory of critical 
phenomena to a quantum phase transition. Hertz ||^ derived an order parameter field 
theory for the ferromagnetic quantum phase transition by considering a simple model 
of itinerant electrons that interact only via the exchange interaction in the particle- 
hole spin-triplet channel. Hertz analyzed this order parameter field theory by means of 
renormalization group (RG) methods. He found a continuous phase transition whose 
critical behavior in the physical dimensions d = 3 and c? = 2 is mean field-like, since 
the dynamical critical exponent z = 3 decreases the upper critical dimension from 
d+ = 4 for the classical case to (i+ = 1 in the quantum case. Despite the somewhat 
artificial nature of this model, it was believed for a long time that the qualitative 
features of Hertz's analysis, in particular the fact that there is mean field-like critical 
behavior for all d > 1, apply to real itinerant quantum ferromagnets as well. 

Here we will show, however, that this belief is mistaken. The properties of the 
ferromagnetic quantum phase transition are much more complicated since the magne- 
tization couples to additional, non-critical soft modes in the electronic system. Math- 
ematically, this renders the conventional Landau- Ginzburg- Wilson approach invalid 
since an expansion of the free energy in powers of the order parameter does not ex- 
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ist. Physically, the additional soft modes lead to an effective long-range interaction 
between the order parameter fluctuations. This long-range interaction, in turn, can 
change the character of the transition from a continuous transition with mean-field 
exponents to either a continuous transition with non-trivial (non-mean field) critical 
behavior or even to a first order transition like in MnSi. 

The derivation of our theory Q follows Hertz Q in spirit, but the technical details 
are considerably different. We consider a d-dimensional continuum model of interact- 
ing electrons, and pay particular attention to the particle-hole spin-triplet contribution 
Q to the interaction term in the action, Sf^^, whose (repulsive) coupling constant we 
denote by Tt. Writing only the latter explicitly, and denoting the spin density by rij, 
the action reads, 

S^Sq + = So + iTt/2) J dx n,{x) ■ n,{x) , (1) 

where Sq contains all contributions to the action other than Sf^^^. In particular, it 
contains the particle- hole spin-singlet and particle-particle interactions, which will be 
important for what follows.^ J dx — J dx Jq^"^ dr, and we use a 4-vector notation x = 
(x, r), with X a vector in real space, and r imaginary time. Following Hertz, we perform 
a Hubbard-Stratonovich decoupling of Sf^^ by introducing a classical vector field M(a;) 
with components that couples to ns{x) and whose average is proportional to the 
magnetization, and we integrate out all fermionic degrees of freedom. We obtain the 
partition function Z in the form 

Z = e-^"/'^ J D[M] cxp[-^M]] , (2) 

where Fq is the non-critical part of the free energy. The Landau-Ginzburg- Wilson 
(LGW) functional $ reads 

$[M] = i y i- S{x - y) M{x) ■ M{y) (3) 

+ ^ On I dxi ... dxn x!".^..i„(a;i, . . . ,Xn)M'-^{xi) . . . M'"(x„) , 

n=2 

where a„ = (— l)"+^/?i!. The coefficients X^") in (I) are connected n-point spin density 
correlation functions of a reference system with action S'o ||. The particle-hole 
spin-triplet interaction Tt is missing in the bare reference system, but a nonzero Tt 
is generated perturbatively by the particle-particle interaction contained in Sq. The 
reference system then has all of the characteristics of the full action 5, except that it 
must not undergo a phase transition lest the separation of modes that is implicit in 
our singling out Sf^^ for the decoupling procedure breaks down. 

is the spin susceptibility of the reference system. Performing a Fourier trans- 
form from X — (x, r) to g = (q, ^2) with wave vector q and Matsubara frequency 11, 
we have for small q and 17, 

x(^)(q ,»)^Xo(q)[l-|»l/|q|] , (4) 

^We note, however, that even in Hertz's original model, where So describes free electrons, interac- 
tion terms get generated upon renormalization. The traditional mean-field results therefore are not 
correct for this model either. 
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where q and are being measured in suitable units, and Xo(q) is the static spin 
susceptibihty of the reference system. If we take the susceptibiHty to be that of a 
non-interacting Fermi gas, xo(q 0) — cq — C2q^, we obtain Hertz's theory. However, 
in a real Fermi liquid, the static spin susceptibility at T = is a non-analytic function 
of q since the magnetization modes couple to additional soft modes, viz. particlc-holc 
excitations in the spin-triplet channel with a ballistic dispersion relation |^ . This non- 
analyticity is crucial for the non-trivial physics of the itinerant ferromagnets discussed 
in this paper. For small wave vectors the static susceptibility is of the form 

Xo(q^0) = co-Cd-i|q|'^-i-C2q2 . (5) 

Here cq, Cd-i and C2 are constants. This holds for l<(i<3;ind = 3 the non- 
analyticity is of the form — C2q^ In |l/q|. Note that all these singularities only exist 
at zero temperature and in zero magnetic field since both a finite temperature or a 
magnetic field gives the particle-hole excitations a mass. 

Using and with = J2q'^J2in^ ^he Gaussian part of <f> can be written, 

$(2)[M] = /"M(g)[io + Cd_i|q|'*-i-f C2q' + cn|f^|/|q|]M(-g) . (6) 

Jq 

Here to = 1 ~ ^tX^^H'i ^ 0,ujn — 0) is the bare distance from the critical point, and 
Co is another constant. 

For the same physical reasons for which the non-analyticity occurs in (^, the 
higher coefficients x''"'* {'n > 2) in (||) are in general not finite at zero frequencies 
and wave numbers. Generally, the coefficient of |]V[|" in $ for |p| — > behaves like 
= ti(") IpI'^+i-". This implies that <i> contains a non-analyticity which in our 
expansion takes the form of a power series in |Mp/|pp. Consequently, the free energy 
functional (||) is mathematically ill defined. However, we will nonetheless be able to 
extract a considerable amount of information. 

The sign of the non-analyticity in the Gaussian term merits some attention since 
it will be responsible for the qualitative features of the ferromagnetic quantum phase 
transition. Perturbation theory to second order in Ft yields Cd-i < ^. This is the 
generic case, and it is consistent with the well-known notion that correlation effects in 
general decrease the effective Stoner coupling However, Ref. has given some 
possible mechanisms for Cd-i to be positive at least in some materials. 



3 Phase transition scenarios 

Depending on the sign of the non-analyticity in the Gaussian term (^) of the free 
energy functional the properties of the ferromagnetic quantum phase transition will 
be qualitatively different. 

We first discuss the generic case of Cd-i < 0. Here the free energy reduces with 
increasing q from zero which implies that a continuous transition to a ferromagnetic 
state is impossible at zero temperature. Two possible scenarios for the phase transi- 
tion arise for Cd-i < 0. The first scenario is based on the observation that a finite 
thermodynamic magnetization m = (|]V[(a;)|), which acts similarly to a magnetic field. 
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cuts off the singularities in the coefficients of the order parameter field theory. There- 
fore, the non-analyticity in x^"^^ leads to an analogous non-analyticity in the magnetic 
equation of state, which takes the form 

tm - Vd'm'^ + um^ = H (d < 3) , (7) 
im - W3m^ ln(l/m) + UTO^ = H (d = 3) , (8) 

where t tunes the transition and u, Vd and are positive constants. H denotes the 
external magnetic field. This equation of state describes a first-order phase transition 
since the next-to-leading term for small m has a negative sign. We have investigated 
this scenario in some detail 0. Since the non-analyticities in x^'^^ and the equation of 
state are cut off by a finite temperature, the transition will be of first order at very 
low T but turn second order at higher temperatures. The two regimes are separated 
by a tricritical point. 

The second possible scenario for the quantum phase transition arising if Cd-i < 
is that the ground state of the system will not be ferromagnetic but instead a spin- 
density wave at finite q. This scenario has not been studied in much detail so far, 
but work is in progress. It is tempting to interpret the spiral ordering in MnSi as a 
signature of this finite-g instability. This is, however, not very likely since a finite-g 
instability caused by our long-range interaction will be strongly temperature dependent 
due to the temperature cutoff of the singularities. As mentioned above, experimentally 
the ordering wave vector is essentially temperature independent. Further work will be 
necessary to decide which of the two possible scenarios, viz. a first-order ferromagnetic 
transition or a continuous transition to modulated magnetic order, is realized under 
what conditions. Moreover, let us point out, that in d = 3 the non-analyticity is only 
a logarithmic correction and would hence manifest itself only as a phase transition at 
exponentially small temperatures, and exponentially large length scales. Thus, it may 
well be unobservable experimentally for some materials. 

We now turn to the second case, Cd~i > which can happen, if one of the conditions 
discussed in Ref. |^ is fulfilled. In this case the self-generated long-range interaction 
is ferromagnetic. Consequently, the ferromagnetic quantum phase transition will be 
a conventional second order phase transition, which can be analyzed by standard 
renormalization group methods. A tree level analysis shows that the Gaussian theory 
is sufficient for dimensions d > = I since all higher order terms are irrelevant. We 
are therefore able to obtain the critical behavior exactly, yet due to the long-range 
interaction it is not mean field-like. The results of this analysis Q can be summarized 
as follows. The equation of state close to the critical point reads 

tm + Vdm'^ -f um^ ^ H (rf < 3) , (9) 
tm ^ vj,m^ / m) + um^ = H (rf = 3) , (10) 

Again, u and v are positive constants. Note the different sign of the non-analytic term 
compared to (||). From (^,^0|) one obtains the critical exponents (3 and 5, defined by 
m ^ t^ and m ^ H^^^ , respectively, at T = 0. For f3 and 6, and for the correlation 
length exponent i^, the order parameter susceptibility exponent 77, and the dynamical 
exponent z, we find 



I3^u = l/{d-l), r] = 3-d, 5 = z = d , (l<d<3) 



(11) 
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and (3 = V ^ 1/2, 77 = 0, 5 = z = 3 for d > 3. These exponents 'lock into' mean-field 
values at d = 3, but have nontrivial values for d < 3. In d = 3, there are logarithmic 
corrections to power-law scaling. Eqs. (|g|,p^) apply to T = 0. At finite temperature, 
we find homogeneity laws for m, and for the magnetic susceptibility, Xm, 

m{t,T,H)^b-^'''m{th^'\Th'*''\Hb^^''') , (12) 

Xrn{t, T, H) = b^/^^Xmitb^^", Tb'^^", Hb'^''') , (13) 

where b is an arbitrary scale factor. The exponent 7, defined by Xm ~ &i T — 
H = and the crossover exponent (j) that describes the crossover from the quantum 
to the classical Heisenberg fixed point (FP) are given by 

7 = /3(^-l)-l , (f> = i^ , (14) 

for all d > 1. Notice that the temperature dependence of the magnetization is not 
given by the dynamical exponent. However, z controls the temperature dependence 
of the specific heat coefficient, — cy/T, which has a scale dimension of zero for all 
d, and logarithmic corrections to scaling for all d < 3 ||ll| , 

-iv{t,T,H) = e{-i~d)\nb + -iv{tb^''',Tb',Hb^^''') . (15) 

Eqs. (^) - (^5|) represent the exact critical behavior of itinerant quantum Heisenberg 
ferromagnets for all d > 1 with the exception of d = 3, where additional logarithmic 
corrections to scaling appear 



4 Influence of disorder 

In this section we briefly discuss the influence of quenched non-magnetic disorder on the 
ferromagnetic quantum phase transition. An approach along the lines of the one for the 
clean case sketched in Sec. || has been developed in Ref. and the resulting effective 
theory is very similar. Again, the magnetization couples to additional soft modes 
(here with diffusive dynamics) which leads to an effective long-range interaction. The 
singularities are even stronger than in the clean case, but they have the opposite sign 
so that the long-range interaction is generically ferromagnetic. Thus in the presence 
of disorder there will be a competition between the ballistic and diffusive singularities, 
and the temperature which cuts off both. For weak disorder the first-order transition 
will survive, while larger disorder leads to a continuous transition. As shown in Ref. 

, the phase diagram becomes very rich, showing several multicritical points and even 
regions with metamagnetic behavior (see Fig. |^). The properties of the continuous 
quantum phase transition occurring for stronger disorder can again be analyzed by 
standard renormalization group methods. It turns out that as in the clean case the 
Gaussian theory is sufficient since all higher order terms are irrelevant. The resulting 
critical exponents are 

7 = 1 , (16) 
for all d > 2, 



j/=l/(d-2) , 77 = 4-d , z = d 



(17) 



600 



Ann. Phys. (Leipzig) 8 (1999) 7-9 





G=Gtce=0.0173| 


G=0.07| 


G=GcE= 


0.092I 


G=0.13<Gc| 


G=0| ^TCP^^ 


^TCEP^ 


CP? 

CEP^ 


\, . 


CF 


CEP 


A" 





-0.02 0.02 -0.02 0.02 -0.02 0.02 -0.02 0.02 -0.02 0.02 
t t t t t 



Fig. 2 Phase diagrams of disordered itinerant ferromagnets in the T-t-plane showing first 
order (solid) and second order (dashed) transitions. G is a dimensionless measure of the 
disorder, and CP, CEP, TCP, and TCEP denote critical points, critical endpoints, tricritical 
points and tricritical endpoints, respectively. See Ref. M| for more details. 



for 2 < d < 4, while = 1/2, 77 = 0, and z = 4 for d > 4. In addition to (i = 4, d = 6 
also plays the role of an upper critical dimension, and one has 

^ = 2/(d"2) , 5^d/2 , (18) 

for 2 < d < 6, while (3 = 1/2, (5 = 3 for d > 6. 

An important problem in disordered systems that has attracted a lot of attention 
within the last years is the influence of rare disorder fluctuations on the properties of 
the phase transition. In the conventional perturbative approach |^ the rare regions 
are neglected. We have developed a theory which includes the rare regions into a 
renormalization group approach. This theory is discussed in a separate paper in this 



volume 1 13 . Here we quote only the final result for the itinerant ferromagnet: Due 
to the effective long-range interaction which stabilizes the Gaussian theory, rare re- 
gions do not change the properties of the ferromagnetic quantum phase transition (in 
contrast to, e.g., itinerant antiferromagnets) . 



5 Conclusions 

To summarize, we have discussed the ferromagnetic quantum phase transition in itin- 
erant electron systems. It has been shown that the critical magnetization modes 
couple to additional, non-critical soft modes in the electronic systems which results in 
an effective long-range interaction between the magnetization fluctuations. We have 
discussed several possible scenarios for the ferromagnetic quantum phase transition. 
In clean systems the generic scenario is a first-order phase transition. This provides 
us with a complete explanation for the nature of the transitions observed in MnSi, 
which in Ref. |^ were attributed to a band structure feature characteristic of MnSi. 
While this feature may well be sufficient to make the transition in MnSi of first order, 
the present theory leads to the surprising prediction that the first order transition is 
generic^ and thus should be present in other weak clean itinerant ferromagnets as well. 
Our theory further predicts in detail how the first order transition will be suppressed 
by quenched disorder. Observations of such a suppression, or lack thereof, would be 
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very interesting for corroborating or refuting the theory. Semi-quantitatively, the the- 
ory predicts that the T region that shows a first order transition will be largest for 
strongly correlated systems. Conversely, since the dependence of the tricritical tem- 
perature on the system parameters is exponential, in some, or even many, systems the 
first order transition may take place only at very low temperatures. This may explain 
why in ZrZn2 no first order transition has been observed H, although the experiment 
does not seem to rule out a weakly first order transition jl2] . 

Under certain conditions, the transition in clean itinerant ferromagnets can also be 
of second order. These conditions have been formulated mathematically. Their precise 
physical or experimental nature has not been investigated so far, but it is known that 
strong correlations are a necessary condition. If these conditions are fulfilled, then 
the critical behavior is known exactly, albeit it is not mean field-like. If the system 
contains sufficiently strong quenched disorder, the transition is always of second order, 
and the non-mean field-like critical behavior has been determined exactly. 

Finally, it has been shown that the presence of rare regions or local moments in the 
disordered case does not change these results, in contrast to, e.g., the case of itinerant 
antiferromagnets. The reason lies in the effective long-range interaction between the 
order parameter fluctuations in the ferromagnetic case, which is sufficient to suppress 
all fluctuations at the critical point, including the static disorder fluctuations respon- 
sible for rare regions. The same suppression of fluctuations is also the reason behind 
our ability to determine the critical behavior for bulk systems, and indeed for any 
dimension, exactly. 

We thank G. Lonzarich and C. Pfleiderer for helpful discussions, and the Aspen Center for 
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